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ABSTRACT 

By means of an automatic digital computer solutions were ob~ 
tained for a number of second order servo problems with backlash in the 
gear train. A complete range of solutions were obtained with variations 
in load friction, motor friction, load inertia, motor inertia, damping co- 
efficient for the combined operation, and backlash angle. These data 
are presented in chart form suitable for the servo design engineer to 
estimate performance of an arbitrary second order system with backlash. 

The program was conducted on the Control Data Corporation Model 
1604 Computer, and this program is presented with explanations of the 
manner in which the differential equations were solved. From data 
generated by the computer phase plane plots are drawn as a further 
illustration of the results of this thesis, and a sample problem is in- 
cluded to show application of the design charts. 

The writer wishes to express his appreciation to Dr. George J. 
Thaler of the Electrical Engineering Department and to Dr. E. J. Stewart 
and Mr. E. N. Ward of the U. S. Naval Postgraduate Computer Center 


for their assistance and encouragement in completing this work. 
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ae Introduction 

In the engineering analysis and application of servomechanisms 
(hereafter called servos) which include gear coupling, the problem of 
backlash must be considered and is often a source of trouble. In theory 
it is desirable to know the system performance of a servo with a given 
distribution of inertia and friction and to accurately determine the size 
of the limit cycle for any given set of variables. In practice, knowledge 
of the performance of a system may be desirable prior to the construction 
of a mock-up or simulation system. 

Existing methods forthe analytical solution of the generalized 
backlash problem are laborious and complicated and lack of the straight- 
forward approach needed in the preliminary design of a servo. In this 
thesis design charts are presented for the second order servo with back- 
lash to accommodate an arbitrary set of system variables. 

Previous work in this problem has been done at the U. S. Naval 
Postgraduate School by Lutkenhouse in reference a wherein the backlash 
problem was solved by construction of the phase trajectory in the graphical 
phase plane. This thesis extends the work of Lutkenhouse. It is recog- 
nized that the construction of a phase trajectory for a given problem 
requires many hours of tedious calculation and graphical construction. 
Thus, in the collection of a large volume of data required for the con- 
struction of design charts for the second order servo with backlash an 


automatic digital computer was employed. 





The work of this thesis first required derivation of the system 
equations, then a computer program was developed to adequately des- 
cribe the equations. The burden of the work subsequently fell in pro~- 
cessing sufficient data to construct design charts to satisfy the require- 


ments of the problem. 





a Equations for the Servo Motion 

In the construction of a set of design charts it is usually desirable 
to provide sufficient flexibility so that the charts may be used for any 
given servo. Thus, equations are presented in Appendix A for the de- 
termination of the system parameters necessary to enter the design charts 
of this thesis for three types of second order servos: the armature control 
motor, the field control motor, and the two phase motor. 

All of these cases behave in a similar manner when coupled to a 
load through a gearbox which has backlash. Initially there is a period 
when the backlash is taker up, and the system performs as a combined 
system. The differential equation describing this motion for an undamped 
natural frequency of one is 

G+ BO +O=1 
where B is equal to twice the damping coefficient, ¢ 

At such point in the cycle of motion wherein the load velocity 
exceeds that of the motor, backlash starts, and the motor is not directly 
coupled to the load. There are two differential equations describing 
this motion. The equation for the coasting load is 

4. +DG.=0 
where D = {e 
The equation for the motor is 


EGy FF by FE. =1 


and F=F... 


where E = ua M 





This motion in the backlash region continues until the total amount of 
backlash is taken up, and at that point gear contact is made. The action 
of the servo upon contact is not a simple one. This motion can follow one 
of several types of contact where the basic laws governing the motion are 
the law of conservation of momentum, and the law of conservation of 
energy which must both be considered. In this thesis it is assumed that 
the contact is a plastic one with no gear bounce; therefore, there is energy 
loss in the gearbox due to the plastic contact and the law of conservation 
of momentum was applied. Following the gear contact phase the system 
is again combined, and the cycle of events beginsagain. 

The characteristics of the second order servo with backlash is 
independent of the undamped natural frequency as shown by the develop- 
ment in Appendix A. Therefore, all of the problems worked in this thesis 
and used in the construction of design charts are based on an undamped 
natural frequency of one radian per second, and the step input is equal 
to one radian. Units of radians and seconds are used throughout. 

The amount of backlash used in the program was 0.3 radian, and 
it was a purpose of the thesis to determine the effect of other magnitudes 


of backlash. 





Si. Solution of the Equations 

The differential equations for the second order servo with backlash 
were solved on the Control Data Corporation Model 1604 Computer. This 
computer is described in reference b, but briefly, it is an all-transistorized , 
stored program, general-purpose digital computer having a large storage 
capacity and exceedingly fast computation and transfer speeds. Input to 
the computer was by means of magnetic tape and punched paper tape, and 
the output data was recorded on magnetic tape. The data was then printed 
on the International Business Machine 717 high speed printer. 

Figure 1 shows a schematic of the problem. In programming for the 
solution of the equations it was necessary to successively compute a 
point for the combined motion, test this point to determine if the load 
would drift away from the motor, then if so, continue the cycle into the 
backlash region. A second test was then introduced to determine the end 
of the backlash region, and if the backlash was fully taken up the contact 
part of the program was initiated. This complete program was accomplished 
by the program of Table I. 

Several instructions in the program are worthy of mention. First, 
it is noted that a subroutine was used in the actual solution of the dif- 
ferential equations. This subroutine based on the Runge-Kutta Gill numer- 
ical integration is contained in a library tape in the U. S. Naval Post- 
graduate Computer Center, and programming instructions are contained 
in reference c. The program was assembled in accordance with the "AR" 
assembly routine of references d and e which also contain a full des- 


cription of each instruction used in this program. 
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In the Runge-Kutta Gill numerical integration four successive 
approximations are used in the computation of a point (in this thesis 
points are 0.01 second apart). However, it was not desired to print out 
each point, so a count routine was inserted in the program to print each 
eighth point. Thus in the actual printout, an example of which is con- 
tained in Table II, test points are at intervals of .08 seconds. The 
printout format shows the number of the point in column 1, the time in 
colian 2, the motor speed in column 3, the motor position in column 4, 
the load speed in column 5, and the load position in column 6. The 
number following the figure printed represents 4 decimal exponent to 
locate the decimal point. 

The remaining item of the program requiring explanation is the 
portion ot the program named switch. This routine is inserted to accom- 
modate the various changes in sign caused by the backlash non-linearity . 
A close study of this sub-routine will show that it was not only necessary 
to change the sign of some of the instructions, but it was also necessary 
to change some of the jump instructions caused by multiple sign changes 
in the backlash test. 

With the basic program of Table I stored in the computer, it is 
merely necessary to insert the desired syster parameters for the parti- 
cular problem to be solved into the computer locations beginning at 20156 


in order to obtain the solution for any desired servo problem with backlash. 





A total of 120 runs for the second order servo with backlash were 
inade and the limit cycle for these problems were determined. Each ruin 
reguiréd approximately three minutes on the computer and ylelded one 
test point suitable for presenting in this report. These test points are 
plotted in Pigures 2 through 11. Figures 12 through 15 are phase trajec- 
tories tor a number of selected runs plotted from an actual printout. 

Figure 2 is a plot of limit cycle versus friction ratio for a damping 
coefficient of 0.1 and inertia ratios of 1/9, 1/4, 1, 4 and 9. 

Figure 3 is a plot of. limit cycle versus friction ratio for the load 
inertia equal to the motor inertia and for damping coefficients of 0.1, 
Omz, 0.3, OF, OWS, and 0.6. 

Figure 4 is a plot of limit cycle versus backlash angle for load 
inertia equal to motor inertia and for a friction ratio of 0.2. This plot 
shows a linear variation for damping coefficients of 0.1, 0.2, 0.3, 0.4, 
0.5, and 0.6, 

Figure 5 is a plot of limit cycle versus backlash angle for load 
inertia equal to motor inertia and for a damping coefficient equal to 0.1. 
This line represents all friction ratios since the Iu - Jy variation of 
Figure 2 is a non-varying straight line. For a damping coefficient other 
than 0.1 the curve would spread to form a pencil of lines. 

Figure 6 is a plot of limit cycle versus backlash angle fora 
damping coefficient of 0.1 and for a friction ratio of 0.2. Variations 


of the curves of this show inertia ratios of 1/4, 1 and 4. 





Figures 7 through 11 are plots of limit cycle versus inertia ratios 
for friction ratios of 0.1, 0.2, 0.4, 0.6, and 1.0. Each plot shows 
curves for variations of damping coefficient of 0.1 through _ 5 

Figures 12 through 15 are actual plots of the phase trajectories as 
plotted from the print-out sheets. From these plots can be determined the 
output Speed and position with the corresponding time. 

Figure 12 shows a case with friction ratio of 0.2 at an inertia 
ratio of .25 and a damping coefficient of 0.1. Figure 13 is at the same 
friction ratio of 0.2 but shows an inertia ratio of 4 and a damping co- 
efficient of 0.6. Figure l4*and 15 illustrates the effect of increasing 
the damping coefficient from 0.1 to 0.2 with load inertia equal to motor 


inertia and at a friction ratio of 1.0. 





Ore Discussion 

After the accumulation of the data it was desirable to present the 
information in the most useful and practical form. The size of the limit 
cycle as shown in the figures of this thesis was previously defined as 
deviation from the equilibrium or average value rather than peak to peak 
variation. The data were plotted in various forms in an attempt to find 
the most useful form. Figure 2 is one such presentation where the data 
is plotted for constant damping coefficient versus friction ratio. Figure 
2 was not found to be the optimum form because in subsequent plots at 
different values of damping coefficient there was appreciable crossing 
and recrossing of lines of constant inertia ratio, and it was found that 
interpolation between these lines was difficult. Further, there was not 
a wide variation of limit cycle with friction ratio as evidenced by the 
level curves of Figure 2. In fact, for the particular case of balanced 
load and motor inertia and damping coefficient there was no variation 
in limit cycle with friction ratio. This plot is shown as a logarithmic 
plot because it was found that the digital program required a number of 
modifications to satisfy the zero friction case. A single case was run with 
zero friction ratio, and the limit cycle was found to be substantially equal 
to that for a friction ratio of 0.04 which is used in this thesis as repre- 
sentative of the zero friction ratio case. 

Figure 3 is a composite plot showing all values of damping co- 
efficient up to 0.6. Figures 4, 5, and 6 are presented to show that the 


Size of the limit cycle is directly proportional to the amount of backlash. 





several runs were made with changes in only the backlash to conclusively 
prove this characteristic. Figure 4 shows the variation with damping co 
efficient with all other variables constant. Figures 5 and 6 are similar 
plots with variations in other factors. This linear variation of limit cycle 
with backlash is significant since it means that design charts can be 
drawn for any constant amount of backlash and variation from this value 
is in direct proportion. In this thesis all charts are for a constant amount 
of backlash equal to 0.3 except for those points in Figures 4, 5, and 6 
where backlash in a variable. 

Figures 7 through 11 @re the design charts found to be the most 
useful form for design estimations. In these charts the size of the limit 
cycle is plotted versus inertia ratio. The logarithmic scale was chosen 
for two reasons: first, it avoided the zero cases found at the extreme 
ends of the scales and for which a separate program would have been 
required, and secondly, this scale gives a spread to the curves in the 
region of inertia ratio equal to one which appears desirable. 

The design charts presented herein are most useful for the range of 
variables encompassed by the test points. While it is believed that some 
estimates can be made from these charts tor systems with very large and 
very small inertia ratios and beyond the limits of the curves shown, best 
results are obtained for values where interpolation between the lines of 
constant damping coefficient is feasible. Appendix B gives a sample 
problem using the design charts. Appendix C gives instructions in the 


use of the design charts. 
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0. Conclusions 

From the results obtained in the accomplishment of this thesis 
it is concluded that the size of the limit cycle of a second order servo 
with backlash can be represented by design charts, one form of which 


are presented in Figures 7 through 11 of this thesis. 


7. Recommendations 

Inasmuch as it was necessary to complete this thesis during a 
Summer period, time was of essence, and the scope of the program was 
definitely limited. Therefore, the following recommendations are sub= 
mitted as a reasonable on of the program. 

1. Rewrite the digital program to include a test to determine the 
actual size of limit cycle for display on the console with no printout 
required. This will permit longer runs with no time required for printout 
of data. Printout capability should still be retained by means of a Jump 
Key. 

2. Extend the program to include the zero friction and zero inertia 
cases. 

3. Extend the program to include a larger range of inertia ratios. 

4, Reprogram the contact phase to include the case of elastic 


contact and to include the “lossy” case for all types of contact between 


elastic contact and plastic contact. 


Il 
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GU 


BACKLASH 


PACTOR2 


IN DEX 


INCRONE 


BLANK 


TEN 


1 


LDA 
STA 
ENA 
STA 
STA 
OTA 
STA 
SLJ 
ENI 
SLJ 
ENI 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
OCT 
OCT 
OCT 
OCT 
OCT 


DEC 
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0 THeELAL 

0 THETA 

0 0 

0 QU 

0 QTHETA 

0 QW 

0 QTHETAL 

4 PRINT 

0 0 

0 SWITCH 

0 0 
oA B equals twice the 

damping coefficient 
-125 D equals the load 
time constant 

22 E equals motor inertia 
=o F equals motor friction 
we G equals motor inertia 
8 H equals load inertia 
l. Integer 
0 Constant initially zero 
23 Amount of backlash 

0100000000000000 
0 
1 Used in count 
0 
10 Used in count 
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20177 
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20202 


20203 


20204 
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20211 


20212 


20213 


20214 


20215 


20216 


20217 


20220 


20221 


20222 
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25075 
21727 
27334 
72702 
35075 
ZR27 
36314 
46314 
37534 
TOPOZ 
44365 
07534 
45075 
2127 
45605 
34121 
46314 
46314 
47024 
60507 
47334 
72702 
94121 
02436 
54365 
07534 
94631 
14631 
99075 
21727 
95341] 
27024 
99605 
34121 
56050 
41217 
96314 
46314 
97024 
60507 
57534 
72702 
64000 
00000 
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13 


14 


15 


16 


17 


18 


ig 


20 


ae 


24 


ao 


DEC 
DEG 
DEC 
DEC 
DEC 
DEC 
DEG 
DEC 
DEC 
DEG 
DEC 
DEC 
DEG 
DEC 
DECq 
DEG 
DEC 
DEC 
DEC 
DEC 
DEC 
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nO 
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.04 


aus 


.06 


«07 


.08 


09 


0 
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12 
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14 
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Z6 


SY 
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1g 


AL, 


eate 
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20237 
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2024 | 
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20243 


20244 


20245 


20246 


20247 
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02436 
64365 
07534 
64631 
14631 
65075 
21727 
65463 
31463 
65605 
34121 
66314 
46314 
67024 
60507 
67146 
63146 
67534 
(PAT NO 
04000 
00000 
04314 
46314 
04631 
14631 
05146 
63146 
05463 
31463 
06314 
46314 
07146 
63146 
14000 
00000 
14631 
14631 
15463 
31463 
16314 
46314 
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32 


35 


36 


40 


44 


45 


48 


90 


39 


60 


65 


70 


80 


90 


100 


120 


140 
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DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 
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. 28 


. 30 
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FOO 


~ 36 


.40 


44 


»45 


. 48 
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OS 


7 


NOS 


aA 


.80 


aSAU 
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20250 20017146 £42180 DEC BR oN, 
31 4 63146 
2029.1 20 0 24000 200 DEC 208 
00 0 00000 
ZOO 00 0 00000 BUF BSS 6 
00 0 00000 
DECO EQU 71000 
RUNGE EQU 60200 
END 
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iggse009- 


Pov I99§ 
23099999 
31999999 
39999999 
47999999 
pov I999 
5 SUSI be, 
71999999 
LIgg9999 
ois ep pe 
qe JI0I999 
10399999 
HL IS999 
Pideoso 9 
iZ/a0009 
13599899 
14369999 
1 198999 
is739999 
Kez oso9'9 
17589999 
ree 99999 
Eggo 99 
Begg Ja99 
20799999 
Zo 99 
22399099 
2299999 
ZeI9I99 99 
Ze LI9999 
24799999 
25990099 
26399999 
YAY Mes ho be he ke 
2/3009 99 
23799999 
ZISOOI99 
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MOTOR 
SPEED 
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15278417 
22518046 
209387916 
35856424 
41895488 
47480616 
52990944 
37 209252 
61321964 
64919110 
67994283 
70544563 
72570424 
(407-5629 
75067101 
79994775 
75551450 
75072608 
73561932 
70007602 
64634564 
57670357 2 
49317406 
39753812 
29138106 
17610490 
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48719261 
44993108 
39850916 
34619795 
29335478 
24032905 
18746033 
13507650 
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12247439- 
27369046- 
48177680- 
74303153- 


10543340 
14121485 
18127574 
2252292 
ZLELOROOS 
2320675 
37640704 
43185760 
48913847 
94783141 
60752245 
66780440 
72827914 
78855989 
84814975 
90570699 
95967717 
10086982 
L0SsTs7938 
10872835 
11149055 
11336614 
11428723 
11419529 
11304005 
14276023 
14557350 
14896797 
15194727 
15450572 
15664046 
15835140 
15964112 
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LOAD 
SPEED 
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14489345 
21757968 
28660703 
35165619 
41244294 
46871885 
52027169 
56692568 
60854140 
64501553 
67628038 
70230320 
72308523 
73866071 
74909556 
75448598 
75495690 
75066030 
74326600 
73587038 
72854835 
72129917 
7 141971068 
70701649 
69998156 
69301663 
68612100 
67929398 
67253489 
48719261 
45474833 
40349320 
35131214 
29856282 
24559521 
19274963 
14035491 
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27260698 
323 19002 
376392387 
43184500 
48912797 
54782300 
60751613 
66780013 
72827688 
7885590" 
84831950 
90748447 
96606073 
10240541 
10814705 
LiSSist55 
11945950 
12503145 
13054795 
13600957 
14141684 
14276023 
14557542 
14896996 
19194332 
15450780 
15664257 
15835352 
15964323 
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APPENDIX A 
DEVELOPMENT OF CHART PARAMETERS 
The purpose of this appendix is to present a procedure whereby the 

factors required to enter the design charts of this thesis may be computed 
for a given servo system. Equations for the field control motor, the arma- 
ture control motor, and the two phase motor are developed. The application 
of a given system represented by a differential equation of the form 

4+ 2binA tah =ah 
requires a transformation to the form 

B +26b +0=8 
since the design charts are based on an undamped natural frequency of one. 
This can be accomplished by a time scaling of the problem. Thus, in the 


first equation above if the independent variables are related by tpt =f ¢ 


where - JB AG dt — ° 
P= He ~ ot! dt ~ 

and ign bu 
When these substitutions are made there results the desired form suitable 
for use with the design charts. From this, it is evident that the damping 
coefficient, G , remains unaltered in the transformation, and the design 
charts, while based on an undamped natural frequency of one, are applicable 
to other frequencies as well. 

The second order servo can be represented by a block diagram as 


follows. 
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LACK LAS 
Three types of motors are considered, the armature control motor, the 
field control motor, and the two phase motor. For the period of time when 
the gears are in contact the system is a simple second order system. The 
differential equations for this combined motion can be derived by equating 


the load torque required to the motor torque developed. For the armature 


control motor this results in the following equation. 


K 


LE ge me 
g eae Bra’ ie 


For the field control motor if the field time constant is negligible the 


equation is similar. 


P* Ge? gr’ Tat AO 


In like manner the differential equation for a two phase motor is 


Fr ae Sr KO KR 


2) aa 7 Kes Pe ae 


o At +2. Lp FIL. jy FSL. a 
where 
Ky = motor torque constant 
Ko = motor generator constant 
Ke. = error measurement constant 
R = armature resistance 
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gear constant 


inertia of motor referred to load shaft 


inertia of motor 


inertia of load 


friction constant of motor 


= friction constant of motor referred to load shaft 


= friction constant of load 


Il 


friction of load plus motor friction referred to load shaft 


Ki K3 


(ee 


Ay Ke 


PAR 


Bee 
Or) 


(from characteristic curves of two phase motor) 


oT 
DE; 


(from characteristic curves of two phase motor) 
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Procedure for Armature Control Motor 


KY KE 
PR 
2. Find Ie referred to load shaft. 


3. Find ty = [ 
Vt Ta 


4. Pind Fp from fe Sot 


1. FindK from XK = 


2 
Seeeeind F, from - _ ae 
6. Findé _ Bek) =) 
24 
rind J, /Jgy 
8. Find F./F eee 
L Total Ce as 


9. Use 6, 7, and 8 to select proper design chart and find size of limit 
cycle for 0.3 rad backlash. 


10. Find size of limit cycle for desired backlash from 


LAS 
LUMI OCFCLE = (CHAR 7 Linnie Cree) (Ae 


Procedure for Field Control Motor 


1. through 4. Same as Armature Control Motor 
r ya pe ae 
fs 
s. Fink 89 O = (GOV re 


6. Continue as Armature Control Motor 





Procedure for Two Phase Motor : 
1, ied K, from characteristic curves of motor 

Ze ering Iu referred to load chart 

3. Find qd). = eae 


Tit TL 


a 
Find ea 


5. Find RO from characteristic curves of motor. 
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6. Find ¢ o-( SMPs) 


7. Continue as Armature Control Motor, 
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APPENDIX B 
SAMPLE PROBLEM 
A Elinco FD-162 115 volt D. C. Shunt Motor, 1/125 Horsepower 

4,000 Rpm is coupled through a gearbox to a load. The constants of the 
proportional error servo are as follows: 
Armature resistance, R=282 ohms 
Motor Generator constant, K=0 .233 volts/rad/sec 
Torque Constant, K 728.2 ounce - inch/ampere 
Motor Inertia, i 0.6] ounce-inch@ 
Motor Friction, tn 00064 Sin OEE) see 
Error measurement constant, K4=10 volts /rad 
Gearbox factor, p=7 
Backlash 0.1 rad 


=6 
Load Inertia J, =16 4x10 slug ag 


Load Friction FP, =. 002 ounce-inch/rad/sec 


Solution: Using the procedure for the armature control motor in appendix A. 


Fae ELE AOE, = eee 
l. Aa =a EOE be 


te ae, = 
an Ty = .. _ sara OG */O 6 sug Ere 


(2 (NEF (FEE) 


ye KO eee = /4.5 Rodgec 
3k a | 226 x0% | | 
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= fr LE ¢ 0006F% - 
jinn. (61075) 13,700 ‘oe 
se 
et _HKGhke . (22 ee 


6a 
a an CMI Qefey Fa td ne 


6. VEE i = Vike’ SEE HO a 
g= (Z br TZ. ys (nC OY ENE) ME 


i; a= “dAUO?® _ po 


Tf Be ee 





oo ee 
——— — a= .O77 
yor. (VENI IAP VO 
9. From Figure 7 representing a friction ratio of .1, and for an 
inertia ration of 2.0 and a damping coefficient of .19 the 
size of the limit cycle for 0.3 radian backlash is .25 radian. 


10. The size of the limit cycle for 0.1 radian backlash is 


60 =  o8ZKGICaN | 
ay ee 


? 


limiteveleis- =< — ~ 
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APPENDIX C 
USE OF DESIGN CHARTS 
Figures 7 through 11 are presented as representing design charts 
for the second order servo with backlash. These charts can be used to 
rapidly determine the size of the limit cycle by the following procedure. 

1. If the system variables are not known, compute these by 
the methods of Appendix A. 

2. Choose the chart representing the Load Friction/Total Friction 
desired. 

3. Read the Size of fe Limit Cycle resulting from the intersection 
of the abscissa representing inertia ratio and the desired 
damping coefficient, interpolating if necessary. 

4. The value determined represents that Limit Cycle resulting 
from a backlash of 0.3 radians. Fora given backlash the 
limit cycle is directly proportional, thus divide the Limit 


Cycle found in 3 above by 0.3 and multiply by the desired 


backlash. 


45 























